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OF  A  SMALL  CURRENT-CARRYING  LOOP  OVER  A 
HOMOGENEOUS  EARTH 


James  R.  Wait  and  Randolph  H.  Ott 
ABSTRACT 


The  basic  theory  of  airborne  EM  surveying,  in  the  time  domain, 
is  considered.  Rather  than  resorting  to  tedious  double  numerical  inte¬ 
gration,  a  more  direct  approach  is  adopted.  This  method,  valid  in  the 
quasi-static  regime,  is  illustrated  for  a  homogeneous  flat  earth.  The 
results  exhibit  a  number  of  clear-cut  features  that  are  relevant  to  re¬ 
mote  sensing.  For  example,  a  vertical  co-axial  loop  system  has  a 
desirable  transient  response  from  the  standpoint  of  yielding  conduc¬ 
tivity  data  without  requiring  accurate  information  on  the  height  of  the 
device  above  the  ground. 


1.  INTRODUCTION 

Remote  sensing  of  the  earth's  surface  layers  often  utilizes  the  eddy 
currents  induced  by  an  active  electromagnetic  source  such  as  a  magnetic 
dipole  in  an  aircraft.  The  secondary  fields  are  detected  either  at  the  aircraft 
or  in  a  towed  bird.  While  the  main  objective  in  commercial  surveys  is  to 
locate  metallic  mineralization,  the  secondary  field  response  is  a  measure  of 
the  effective  ground  conductivity  of  the  overburden.  This  in  turn  can  often  be 
related  to  the  rock  strength..  Currently,  there  is  a  great  interest  in  utilizing 
transient  source  signals  rather  than  the  more  prosaic  time  harmonic  fields. 

In  this  note,  we  wish  to  present  the  basic  theory  for  the  time  domain 
response  of  a  magnetic  dipole  or  small  loop  located  over  an  idealized  flat 
homogeneous  earth.  The  results  are  basic  to  a  proper  understanding  of  air¬ 
borne  electromagnetic  geophysical  surveys. 

Z.  FORMULATION 

Our  model  is  very  simple;  we  imagine  the  source  to  be  a  small  loop 
(or  vertical  magnetic  dipole)  of  area  dA  located  at  height  h  over  a  non¬ 
magnetic  conducting  half-space  of  conductivity  o  .  We  choose  a  cylindrical 
coordinate  system  (p,  0,  z)  with  the  dipole  at  z  =  h  and  the  earth's  surface 
at  z  =  0.  We  assume  for  sake  of  simplicity  that  all  displacement  currents  may 
be  neglected.  This  is  justified  by  the  fact  that  all  linear  dimensions  of  the 
problem  are  extremely  small  compared  with  the  significant  free  space  wave¬ 
lengths  in  the  problem. 
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The  time  harmonic  solution  for  a  current  I  exp(i  wt)  in  the  loop  is 
well-known.  The  resulting  magnetic  field  components  at  P  are 


h  (i<o)  = 

P  dpdz 


H  (m  =  jiSiial 

Z  3z2 


where  n(icd)  is  the  Hertz  potential  given  by 


I (id?)  dA 

.  C  _L  1 1 

4ir 

— i 

1  * 

■ 

1 

1  K 

I  (iu?)  dA 

o 

/■»*  -X(h+z) 

\  e  J  M  0\  dX 

4  ir 

\  i  J  \  hp)  QA 

o  ( X2  +  a  |a  i  o>) 2  +  X 

l 

i  _ 
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[p2+  (z+h)2]2  ,  u  =47r  x  10  and 

Bessel  function.  If  now  the  current  in  the  loop  is  I  6(t).  where  6 (t)  is  the 

o 

unit  impulse  function,  we  find  without  difficulty  that  the  transient  field  compon¬ 
ents  are  to  be  obtained  from 


where 


h  (t)  =  d2  Il(t)/dpdz 
P 

h  (t)  =  32  n(t)/3  z2 
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where  £  denotes  the  inverse  Laplace  transform  operator. 
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3.  THE  TRANSIENT  SOLUTION 

Using  Pair  No.  543.5  in  Campbell  and  Foster's  [  1949]  tables,  we  find 
that  (7)  is  equivalent  to 


All(t)  = 


I  dA  , 

O _ £_ 

4n  Jv 


r  j  F(— .)  e  dx 


where 


1  *T  r 

F{T)  =  T2  e"  [l  -(irT) 


4  T 


e  erfc 


(T2)] 


and  where  erfc  is  the  complement  of  the  error  function.  By  changing  the  vari¬ 
able  of  integration  in  (8)  to  x  =[t/(o4)  ]  2  X,  we  can  write 


AH(t)  = 


I  dA 
o  Z 


i,  Q  J  [xe  X  -  ff2  x2  erfc(x)  ]  e  Hx  Jq(Px)  dx 


where  H  =(a  |i /t)2  (h+z)  and  P=(ofi/t)2p.  In  order  to  deal  effectively  with 
the  integral  in  (10),  we  write 


'■HXj0<p*>  =1 


(-1)  m 

a  - - —  x 
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re  a  =(-l)m  (c^V^x”1)  e  J  (Px)  1  I  .In  the  special  case  where 
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P  =  0  (i,  e. ,  p  =0),  we  see  that  a  =H  . 

m 

Using  (11),  we  find  that  (10)  becomes 

I  dA  2  v-'  (-1)  a 

An<tl  ■  V -jTTC-r-)  I  sr 
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Now  it  follows  from  the  basic  definition  of  the  Gamma  function  T(n)  of 
argument  n  that 


C  2n-l  -x2 

J  * 


dx  =  r(n)/Z 
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which  has  the  form  of  the  first  inf  O  g  V/f  1  in  M  ,  A  Iso  t  f  Jnm  ftyr\  a  rijj 

Ryzhik  [  1966],  (No.  6281,  pg.  648),  we  have 


erfc(x)  dx  = 


r(g+i) 

2q^ 


which  is  the  form  of  the  second  integral  in  (12).  Thus 
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In  order  to  obtain  explicit  series  representations  for  the  field  components,  we 
need  to  differentiate  (15)  according  to  (4)  and  (5).  Of  special  interest  is  the 
response  h  (t)  in  the  case  p  =0.  Then  we  find  without  difficulty  that 


I  dA 

h  (t)  =  - - 
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where  H  =  (o  U  /t) 2  (h+z)  .  From  (17),  it  is  evident  that  for  large  values  of  t 


I  dA  .  ..  3/2 


w  j#  4, 

which  indicates  that  tue  magnetic  field  response  varies  as  t  for  sufficiently 
large  times.  This  result  is  in  agreement  with  a  similar  recent  analysis  of  the 
problem  by  Ott  [  1971] ,  but  it  disagrees  with  the  analysis  by  Keller,  et  al  [  1970] , 
who  made  unjustified  approximations  in  the  integrand  of  the  time -harmonic 
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integral  representation.  It  is  remarkable  that  (18)  does  not  depend  explicitly 
on  the  distance  between  source  and  observer.  This,  of  course,  would  not  be 
the  case  in  general.  Nevertheless,  it  suggests,  that  for  a  vertical  co-axial  loop 
system,  the  transient  coupling  should  be  relatively  insensitive  to  height  of  the 
device  above  the  ground. 

The  behavior  of  the  transient  response  at  small  times  t  can  be 
obtained  directly  from  (8)  and  (9).  For  example,  if  t  is  sufficiently  small 
F(T)  T2  and  then 


All(t)  ~ 
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*Jn 


XJ  Up)  e'x(h+z)dx 
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(19) 


This  shows  that  £II(t)  and  the  field  components  h  (t)  and  h  (t)  vary  as  t-2 

P  z 

for  sufficiently  small  times. 

In  the  general  case  where  time  t  is  neither  small  nor  large,  it  appears 
necessary  to  evaluate  numerically  the  integral  of  the  type  given  by  (10).  Some 
examples  from  this  approach  are  given  in  Figure  1  where  we  plot  h  (t)  and 

5  ^ 

h  (t)  normalized  to  4 tt  04  p  /l^dA  as  a  function  of  the  time  parameter  T=t/cr4p2 

for  various  values  of  the  ratio,  R  ={h  +  z)/p.  The  ordinates  are  labelled 

H  (T)  and  H  (T)  to  denote  the  normalization, 
z  p 

The  small  time  and  the  large  time  asymptotes  are  also  shown  on  the 
curves.  These  are  discussed  in  the  next  section. 

4.  DISCUSSION  OF  RESULTS  AND  CONCLUDING  REMARKS 
The  first  two  leading  terms  in  the  large  time  response  are  found  to  be 


H  (T) 
z- 


5VFT 


w 


°Cd7F>  — 


(20) 


This  result  is  compared  with  the  numerical  results  in  Figure  1  for  a  height  to 
separation  ratio,  R,  of  1  fz.  The  results  of  (20)  and  the  numerically  computed 
results  show  very  good  agreement  for  large  times. 
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(21) 


This  small  time  response  is  compared  with  the  numerical  results  in 
Figure  1  for  R  =  l/2  .  Again,  good  agreement  between  the  methods  is  obtained. 
An  interesting  feature  of  the  small  time  response  in  (21)  may  be  noted  by  vary¬ 
ing  R  for  fixed  T.  As  we  vary  R,  there  should  be  a  minimum  in  the  small 
time  response  at  R  =  </3/»l/’2  .  Indeed,  as  we  vary  R  in  Figure  1,  we  find  that 
increasing  R  to  3/2  increases  the  response, and  decreasing  R  to  1  also 
increases  the  response. 

As  indicated  in  Figure  1,  we  find  that  the  numerical  results  for  h  (t) 

z 

show  a  reversal  in  polarity,  while  those  for  h  (t)  do  not.  This  change  in 

P 

polarity  does  not  occur  for  h  (t)  for  values  of  R  greater  than  about 

z 

1.30.  Obviously,  from  Figure  l,the  small  and  large  time  responses  are  not 
sufficient  to  completely  describe  the  response  at  intermediate  times.  This 
points  up  the  necessity  of  dealing  with  the  full  integral  representation  of  the 
problem  if  a  proper  understanding  of  the  transient  phenomena  is  to  be  achieved. 
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Fig.  1  Impulse  field  responses  H  (T)  as  a  function  of  no 

time  for  various  R  values.  H  (T)  is  shown  for 

P 
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